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ABSTRACT 

This  report  outlines  our  research  activities  for  the  period  of  the  AFOSR 
grant.  It  Is  not  Intended  to  give  a  detailed  description  of  our  work.  Such 
details  are  described  In  our  publications  and  cited  references.  The  report, 
however,  defines  the  scope  of  our  research  on  the  unsteady  aerodynamics  and 
the  stability  analysis  of  turbomachine  components  and  Its  relevance  to  ongoing 
technological  developments  In  turbomachine  design. 

The  main  topic  of  our  research  Is  the  unsteady  aerodynamics  of  lifting 
airfoils  subject  to  three-dimensional  gusts.  However,  the  mathematical 
methods  we  have  developed  for  this  purpose,  evolved  into  a  general  theory  for 
unsteady  vortical  and  entroplc  distortions  of  potential  flows.  This  theory 
can  be  also  applied  with In  the  approximation  of  the  rapid  distortion  theory  of 
turbulence  to  study  the  change  In  turbulence  characteristics  near  the 
stagnation  point  of  a  bluff  body. 

The  mathematical  theory  was  applied  to  study  the  aerodynamics  of  lifting 
airfoils  subject  to  three-dimensional  gusts.  For  low  Mach  numbers,  an 
Integral  equation  was  derived  and  a  numerical  scheme  was  developed  for  Its 
solution.  All  previous  gust  analyses  correspond  to  particular  cases  of  our 
present  work.  We  have  therefore  compared  our  results  with  previous  linear  and 
nonlinear  theories.  These  results  show  significant  three-dimensional  and 
nonlinear  effects. 

Finally,  we  have  carried  out  a  comparison  between  our  previously 
developed  theory  for  highly  loaded  cascades  with  data  from  experiments  being 
conducted  by  Frank  Carta  at  United  Technologies  Research  Center.  The  results 
show  excellent  agreement  between  our  analysis  and  the  experimental  data. 
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I.  INTRODUCTION 


The  present  research  Is  mainly  motivated  by  the  ongoing  technological 

developments  in  turbomachine  blading  systems.  The  new  design  trends  in  engine 

technology  tend  toward  increased  thrust  per  unit  engine  weight,  more  fuel 

efficiency,  and  more  engine  stability  at  various  operating  conditions.  The 

performance  of  the  engines  ,1s  greatly  enhanced  by  higher  blade  loading  and 

increased  flow  speed.  This  creates  many  new  problems  in  fluid  mechanics, 

/  _  , 
aerodynamics,  structural  design,  and  aeroelasticity. 

Our  research  relates  to  the  stability  of  the  engine  components  such  as 
blades  and  guide  vanes.  These  structural  components  are  subject  to  heavy 
aerodynamic  loads  and  produce  large  air  turning  inside  the  engine  from  one 
stage  row  to  another.  In  addition  to  those  desired  high  aerodynamic  loads 
acting  upon  rotor  and  stator  blades,  there  are  unwanted  fluctuating  loads  that 
produce  flutter  and  forced  vibrations  of  turbomachinery  blades. 

Indeed,  flows  in  turbomachines  are  highly  nonuniform.  Their  irregular 
patterns  initially  result  from  inlet  distortion  and  inlet  turbulence  as  they 
enter  the  engine.  Subsequently  flow  Interaction  with  rotors  and  stators 
produces  more  irregular  flow  patterns  characteristic  of  each  stage  of  the 
machine.  For  example.  Figure  1  shows  a  schematic  representation  of  flow 
behind  a  rotor  blade  row.  We  notice  the  formation  of  secondary  flows 
migrating  outwardly  from  the  hub,  the  development  of  blade  surface  viscous 
boundary  layers  which  extends  behind  each  blade,  and  the  formation  of  tip  and 
hub  swirling  vortices. 

Figure  2  shows  a  stator  blade  in  the  wake  of  a  rotor.  The  velocity 
profile  behind  the  rotor  has  periodic  Irregularities  due  to  viscous  wakes  and 
swirling  vortices  particularly  at  the  tip.  Because  of  the  periodic  nature  of 
a  rotor  blade  row,  it  Is  proper  to  assume  that  the  velocity  defect  and  the 


swirling  vortices  have  essentially  the  same  periodic  structure  as  the  rotor 
Itself.  This  Irregular  but  periodic  flow  pattern  Is  moving  with  a  velocity  Up 
with  respect  to  a  stator  blade.  The  stator  blades  are  hence  subject  to  a 
vortical  wave  called  a  gust.  This  phenomenon  was  first  noticed  by  Kemp  and 
Sears  [1].  Usually,  the  Irregularities  in  the  flow  are  small  compared  to  the 
mean  flow  velocity.  That  Is,  If  Is  a  characteristic  of  the  velocity  defect 
and  UL  Is  the  mean  flow  velocity,  then  U^/U,,,  «  1.  Thus,  the  gust  Is  a 
vortical  wave  essentially  convected  by  the  mean  flow. 

Figure  3  shows  a  schematic  decomposition  of  the  gust  Into  longitudinal 
(chordwlse)  and  transverse  components.  There  is  also  a  spanwlse  component 
which  is  perpendicular  to  the  plane  of  the  figure.  The  three  components  of 
the  gust  are  all  significant  and  for  a  lifting  airfoil  (loaded  blade)  they  are 
strongly  coupled  together  by  the  mean  flow  around  the  airfoil  as  we  have  shown 
recently  In  [2]. 

The  main  objective  of  the  present  research  is  to  formulate  an  aerodynamic 
theory  for  a  three-dimensional  gust  Interacting  with  loaded  airfoils,  and  to 
develop  mathematical  and  numerical  procedures  to  calculate  the  fluctuating 
pressure  and  forces  acting  upon  those  airfoils. 

Although  our  work  Is  motivated  by  turbomachine  technology,  we  should  note 
that  our  research  Is  of  a  basic  nature  and  aims  at  advancing  and  extending  our 
knowledge  and  understanding  of  unsteady  flows  as  they  interact  with  airfoils 
and  blades. 


II.  GENERAL  OUTLINE  OF  OUR  RESEARCH  ACTIVITIES 


The  primary  objective  of  the  present  research  Is  to  calculate  the 
unsteady  forces  acting  upon  lifting  airfoils  subject  to  three-dimensional 
gusts.  However  as  we  formulated  the  mathematical  problem  we  found  as  a 
spin-off  that  the  theory  can  be  extended  to  a  general  splitting  theorem  for 
the  velocity  field  of  unsteady  vortical  and  entroplc  distortions  of  potential 
flows. 

Hence  we  developed  a  general  mathematical  theory  for  studying  the 
Interaction  of  non-uniform  unsteady  flows  with  bodies.  A  particular 
application  of  the  theory  Is  the  problem  of  a  three-dimensional  gust 
Interacting  with  an  airfoil  or  a  cascade  of  airfoils.  Other  Important 
applications  can  be  found  In  the  field  of  aeroacoustics  and  fan  and  helicopter 
noise.  Our  theory  can  also  be  readily  applied  under  certain  conditions,  to 
Investigate  the  change  In  turbulence  characteristics  particularly  near 
stagnation  points  and  Its  effects  on  boundary  layer  stability. 

We  are  applying  our  mathematical  theory  to  essentially  two  problems. 

First  we  studied  the  three-dimensional  gust  past  an  airfoil  with  camber, 
thickness  and  mean- flow  Incidence.  Second,  we  are  presently  investigating  the 
change  In  turbulence  characteristics  near  the  stagnation  point  of  a  bluff  body 
and  an  airfoil. 

Finally,  we  have  carried  out  a  comparison  of  our  previously  developed 
theory  of  an  oscillating  cascade  of  highly  loaded  blades  with  the  ongoing 
experimental  results  of  Frank  Carta  at  United  Technologies  Research  Center. 
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III.  UNIFORMLY  VALID  SPLITTING  OF  UNSTEADY  VORTICAL 
AND  ENTROPIC  DISTORTIONS  OF  POTENTIAL  FLOWS 

In  many  aerodynamic  problems  one  deals  with  nearly  Inviscld  large 
Reynolds  number  flows  past  solid  bodies.  For  streamlined  bodies  such  as 
airfoils,  the  flow  Is  well  approximated  by  a  potential  flow  except  of  course 
In  the  boundary-layer  and  In  the  limited  region  where  separation  occurs.  For 
bluff  bodies,  the  separated  region  extends  Into  the  wake  and  the  potential 
flow  Is  only  valid  outside  the  wake.  The  problem  we  are  concerned  with  here 
is  when  small  disturbances  are  added  upstream  which  will  produce  distortion  in 
the  flow  field.  Two  Important  cases  arise 

(1)  Large  Structure  Disturbances 

This  Is  the  case  where  the  length  i'  and  the  velocity  u'  defining 
the  scale  of  the  disturbance,  and  the  length  i  and  the  velocity  IL 
characterizing  the  flow  upstream  are  such  that 


I'/z  »1  or  0  (1) 


1 

W 


«  (|j-)  «  1 
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where  Re  *  U»  */v  Is  the  Reynolds  number.  This  is  the  case  of 
disturbances  produced  in  turbomachines  by  rotor-stator  Interaction, 
and  In  general  the  case  referred  to  by  aerodynami cists  as  the  gust 
problem. 

(11)  Small  Structure  Disturbances 
This  Is  the  case  where 
1 

“7  t 1  1  ii 1  » ' 

Re  <<  ^  1  and  ^  U  X  * 
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This  is  the  case  of  small  structure  turbulence  interacting  with  a  body. 
However,  the  integral  scale  is  assumed  larger  than  the  boundary  layer 
thickness  so  that  no  direct  coupling  between  the  incident  disturbance  and  the 
viscous  boundary-! ayer  will  occur. 

In  both  cases  a  linear  inviscid  theory  can  be  developed  to  model  the 

fl  ow. 

If  the  body  is  a  flat  plate  at  zero  incidence,  then  there  is  a  complete 
uncoupling  between  the  mean  flow  and  the  disturbance.  The  mathematical 
treatment  reduces  to  a  Laplace  equation  for  an  incompressible  flow  and  was 
carried  out  first  by  Sears  [3].  For  two-dimensional  compressible  flows  and 
oblique  gusts  the  problem  can  be  reduced  to  solving  a  Helmholtz  equation  as 
shown  by  Graham  [4].  For  a  cascade,  an  integral  equation  formulation  was 
first  developed  by  Lane  and  Friedmann  [5].  Whitehead  [6]  gave  detailed 
numerical  calculations  for  arbitrary  interblade  phase  angle.  Atassi  and  Hamad 
[7]  extended  the  cascade  theory  to  three-dimensional  disturbances  in  subsonic 
f 1 ows . 

For  a  lifting  airfoil  the  coupling  between  the  mean  flow  and  the 
disturbance  is  very  strong  as  shown  by  Goldstein  and  Atassi  [8]  for  a 
two-dimensional  gust  and  by  Atassi  [2]  for  a  three-dimensional  gust. 

For  a  bluff  body.  Hunt  [9]  treated  the  Interaction  of  incident 
disturbance  with  a  circular  cylinder  by  generalizing  the  "rapid  distortion 
theory  of  turbulence"  developed  by  Batchelor  and  Proudman  [10]  and  Ribner  and 
Tucker  [11].  These  theories  predict  changes  occurlng  in  turbulent  flows  that 
are  distorted  by  solid  obstacles. 

Let  us  review  the  mathematical  methods  developed  to  deal  with  this 


general  class  of  unsteady  vortical  disturbances. 
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1.  Traditional  Splitting  of  the  Velocity  Field 

*► 

In  this  case  the  unsteady  velocity  u  Is  split  Into  a  solenoidal 
-► 

component  us  and  an  irrotatlonal  component  7$ 

u  =  u$  +  7*.  (1) 

Introduce  the  vector  &  such  that  7 •%  =0  and  such  that 

u$  =  7xX  .  (2) 

This  leads  to  three  Poisson  equations 

72*  =  -  J  (3) 

where  a>  Is  the  vortlclty  which  can  be  determined  using  Cauchy's  solution 

[12, p.  276]  which  depends  on  the  Lagranglan  coordinates  of  a  fluid  particle. 
The  function  $  satisfies  the  equation 


LU)  s 


°o* 

ITT’ 


i_7 

po 


7-(P  U  ) 


(4) 


The  operator  L  defined  in  (4)  Is  a  second-order  linear  convective  wave 
operator  with  non-constant  coefficients.  Note  that 


Is  the  linearized  material  derivative,  and  p0,  c0  represent  the  density  and 
speed  of  sound  of  the  mean  flow  field,  respectively. 

This  method  leads  to  three  Poisson's  equations  and  an  inhomogeneous  wave 
equation  (4).  For  an  Incompressible  flow  (4)  reduces  to  Laplace  equation. 


2.  Goldstein's  Splitting  of  the  Velocity  Field 


Instead  of  splitting 
Irrotatlonal  fields,  Goldstein 


the  velocity  field  Into  solenoldal  and 
♦ 

[13]  splits  u  Into  a  known  vortical  solution 


(H)  =  Aj  (a,o)  yJ- 


(5) 


and  an  irrotatlonal  field  v<j>  satisfying 


L  (+) 


v - ( pouH ) 


(6) 


Here  a  represents  the  Lagranglan  coordinates  of  a  fluid  particle  and  Aj 
are  the  disturbance  components  upstream. 

The  advantage  of  Goldstein's  splitting  lies  in  the  fact  that  the 
mathematical  problem  has  been  reduced  to  solving  equation  (6).  So  Instead  of 
four  partial  differential  equations  we  have  only  one  single  equation. 

The  method  has  been  applied  successfully  by  Goldstein  and  Durbin  [14]  to 
flows  with  no  stagnation  points. 


3.  Flows  with  Stagnation  Points 

Practically,  all  flows  past  airfoils  and  obstacles  have  a  stagnation 
point  on  the  body.  In  this  case  the  Lagranglan  coordinates  used  in  both  the 
traditional  and  the  Goldstein's  splittings  have  a  singularity.  As  a  result 
the  vortlclty  <d  and  the  vortical  velocity  u^  are  both  singular  along  all 
surfaces  of  the  body  and  Its  wake.  In  fact  we  have  proved  the  following 
result: 

"For  a  flow  past  a  body  with  a  stagnation  point,  both  the  traditional 
and  Goldstein's  splittings  lead  to  a  singular  and  Indeterminate  vortical 
solution  along  the  body  surface.  The"  boundary-value  problem  for  the 
potential  part  of  the  disturbance  would  then  have  to  satisfy  singular  and 
Indeterminate  Neumann's  boundary  conditions  at  the  body  surface." 
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This  singular  and  indeterminate  behavior  all  along  the  surface  of  the 
body  is  a  major  difficulty  for  analytical  and  numerical  solutions  to  the  gust 
problem  and  to  the  analysis  of  turbulent  flows.  In  fact,  we  can  even  state 
that  unless  the  solution  is  regularized,  the  treatment  will  be  incorrect. 

In  order  to  illustrate  the  effect  of  the  stagnation  point,  we  consider  the 
airfoil  shown  in  Fig.  4.  The  mean  flow  streamlines  are  also  shown.  Figures  5 
to  12  show  the  streamwise  and  normal  components  of  the  vorticity  and  the 
vortical  velocity  for  an  initial  harmonic  normalized  disturbance  at  45°  to  the 
mean  flow.  The  reduced  frequency  Iq  *  4.  The  choice  of  a  large  reduced 
frequency  is  made  to  illustrate  the  full  variation  of  the  vorticity  and 
vortical  velocities  over  the  airfoil  surface.  Because  of  the  singular  behavior 
at  the  airfoil  surface,  we  only  considered  the  streamlines  i|»  2  ±  0.01.  The 
abscissa  s  is  the  length  of  the  arc  along  the  airfoil.  The  origin  s  =  0 
correpsonds  to  the  stagnation  point  location,  and  the  other  two  vertical  lines 
correspond  to  the  trailing  edge  on  the  pressure  side  and  the  suction  side.  The 
amplification  of  the  streamwise  component  of  the  vorticity  and  the  normal 
component  of  the  vortical  velocity  are  quite  marked.  They  will  become,  of 
course,  infinite  for  <|>  =  0  as  we  move  closer  to  the  surface  of  the  al  ,oil. 

The  phase  of  all  the  components  will  be  indeterminate  at  the  airfoil  surface. 

Hence,  if  one  uses  the  traditional  method  of  solution,  Eqs.  (3,4)  will 
have  singular  and  indeterminate  inhomogeneous  terms.  On  the  other  hand  if  one 
uses  Goldstein's  method,  the  vortical  solution  uH  Is  singular  and 
Indeterminate  along  the  surface  of  the  airfoil  and  Eq.  (6)  will  have  to 
satisfy  singular  and  indeterminate  boundary  conditions  all  along  the  airfoil 
surface  and  its  wake. 
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4.  Uniformly  Valid  Splitting  of  the  Velocity  Field 

Our  mathematical  procedure  can  essentially  be  summarized  by  the 
following  theorem  : 

"The  velocity  field 


u 


(R) 


:(H) 


v* 


where  $  satisfies  the  equation 


-  0, 


and  the  boundary  condition 

uH  +  jj  *  0,  on  Airfoil  and  Wake, 


Is  regular  and  has  the  remarkable  property  that  its  streamwise  and  normal 
components  vanish  at  the  body  surface  and  Its  wake.  The  singular  and 
Indeterminate  character  of  u^  Is  eliminated." 

Hence  we  propose  the  following  mathematical  method  for  all  problems 
Involving  vortical  and  entroplc  distortions  of  potential  flows. 

u  *  u(R)  +  v$  ,  (7) 


where 


£  U(R)  ♦  U<Ri  .  70 


(8) 
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With  homogeneous  conditions  for  u^R)  along  the  tangent  and  the  principal 
normal  of  the  body  surface  and  its  wake  If  any,  and 

L  U)  -  7-p  u(R)  (9) 

po 

with  the  conditions 

*  0  along  body  surface  and 

4  /  R  j 

v<p  +  u  =  prescribed  upstream  conditions. 

More  details  about  this  splitting  were  given  in  [15]  and  will  be  soon 
published  in  [16]. 

General  expressions  for  u(R)  were  also  derived  for  flows  that  exist  for 
all  times  for  three-dimensional  disturbances  past  two-dimensional  and 
axi symmetric  bodies. 

5.  Incompressible  Flows 

For  an  incompressible  flow  Eq.  (9)  reduces  to 

v2*  =  -  v-u(R)  (10) 

and  hence  the  problem  reduces  simply  to  solving  a  Poisson's  equation. 

For  a  two-dimensional  body.  It  Is  possible  to  consider  a  Fourier 
expansion  spanwise,  and  the  preceding  equation  reduces  to 

(v2  +  k2)  *  »  -  V‘u(R) 


(11) 
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where  v2  is  the  2-D  Laplace  operator  and  IC3  is  the  wave  number  In  the  span 
direction.  This  equation  Is  called  a  Kleln-Gordon  equation  [18]  and  usually 
It  occurs  in  studies  of  oscillations  In  resisting  media.  Thus,  one  should 
expect  a  smaller  fluctuating  pressure  for  three-dimensional  oblique  gust  than 
for  two-dimensional  gusts.  These  results  will  be  further  discussed  in  our 
study  of  the  gust  problem. 

IV.  LIFTING  AIRFOILS  SUBJECT  TO  THREE-DIMENSIONAL  PERIODIC  GUST 


The  gust  problem  is  illustrated 
a  blade  can  be  expanded  in  a  Fourier 
generality,  we  can  treat  the  case  of 

.  t  It«  -  U.t] 
u  *  A  e 


In  Figure  2.  The  disturbance  as  seen  by 
series.  Therefore  without  loss  of 
a  harmonic  disturbance  of  the  form 

(12) 


The  general  theory  can  be  developed  as  an  application  to  our  splitting 
theorem  and  the  mathematical  problem  Is  reduced  to  finding  the  regular 
vortical  velocity  u^,  and  then  solving  equation  (9)  subject  to  upstream 
condition  (12)  and  homogeneous  condition  along  the  airfoil  surface.  In 
addition  one  should  account  for  the  effect  of  vortex  shedding  In  the  wake 
which  here  will  be  represented  by  a  vortex  sheet. 

The  early  treatments  of  the  gust  problem  dealt,  of  course,  with  flat 
plate  airfoils.  These  are  the  so-called  linearized  theories  which  we  briefly 
review  here  since  they  will  be  used  for  comparison  with  our  present  work. 
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1.  Linear  Theories  for  Airfoils  and  Cascades  In  Nonuni  form  Motions 

The  aerodynamics  of  flat  plate  airfoils  subject  to  a  gust  have  been 
extensively  developed  at  all  flow  regimes.  They  constitute  an  Important 
reference  case  for  our  present  work,  and  our  experience  In  linear  unsteady 
aerodynamics  Is  of  great  help  for  pursuing  the  present  nonlinear  problem. 

Although  the  linear  treatment  brings  about  considerable  simplification  of 
the  mathematical  analysis,  the  only  closed  form  analytical  solution  for  the 
fluctuating  lift  and  moment  Is  that  of  Sears  [3]  for  the  flat  plate  in 
Incompressible  flow  subject  to  a  transverse  gust.  For  compressible  subsonic 
flows,  an  Integral  equation  Is  commonly  used  In  the  analysis.  For  a  single 
airfoil  we  have  an  equation  derived  by  Posslo  [18],  and  for  a  cascade  the 
treatment  was  first  made  by  Lane  and  Friedman  [5].  Detailed  numerical 
methods  were  given  by  Whitehead  [6,19],  and  numerical  codes  are  now  commonly 
used  for  two-dimensional  gusts  Interacting  with  flat  plates  and  cascades. 

The  case  of  an  obliquely  propagating  gust  Interacting  with  a  flat  plate 
would  be  the  linear  counterpart  of  the  three-dimensional  gust  Interacting  with 
a  lifting  airfoil  we  propose  to  Investigate.  This  Is  because  the  third 
longitudinal  component  of  the  gust  does  not  induce  any  fluctuating  forces  when 
the  lift  coefficient  of  the  airfoil  Is  zero.  Fllotas  [20],  Mugrldge  [21],  and 
Graham  [22]  using  different  formulations  studied  the  response  function  of  a 
flat  plate  to  an  oblique  gust. 

For  a  cascade  of  flat  plate  airfoils,  Atassl  and  Hamad  [7]  derived 
similarity  rules  for  three-dimensional  gusts  In  compressible  flows.  Atassl 
and  Hamad  also  developed  a  computer  code  for  the  aerodynamics  of  swirling 
flows  Interacting  with  cascades  and  applied  their  results  to  determine  the 
level  of  noise  generated  by  tip  vortices. 


The  most  significant  features  of  linear  theories  are  given  below: 


(1)  Single  Flat  Plate  Airfoil 

The  unsteady  lift  depends  on  the  two  wave  numbers  Iq  In  the 
flow  direction  and  IC3  In  the  span  direction.  Note  that  Iq  is  also  the  reduced 
frequency.  Figure  13  shows  the  real  and  Imaginary  parts  of  the  lift 
coefficient  L  for  the  spanwlse  wave  number  IC3  *  0.0,  0.2,  0.4,  0.6,  1.0,  2.0, 
and  3.0.  The  reduced  frequency  Is  varied  along  each  curve  from  0  to  5.  The 
Mach  number  Is  0.8.  The  wave  number  k2  in  the  transverse  direction  was 
factored  out  since  in  the  linear  case  it  does  not  affect  the  solution. 

Figure  13  underscores  the  Importance  as  well  as  the  complexity  of  the 
three-dimensional  effects.  The  essential  features  of  the  3-D  gust  are: 

•  At  low  reduced  frequency,  the  unsteady  lift  increases 
as  the  mean  flow  Mach  number  increases. 

•  The  obliqueness  of  the  gust  significantly  reduces 
the  magnitude  of  the  unsteady  lift  at  low  reduced 
frequency. 

•  The  lift  function  does  not  decrease  with  increasing 
k3  when  Iq  -  I13. 

(11)  Cascade  of  Flat  Plate  Airfoils 

The  cascade  Introduces  a  large  number  of  parameters  such  as 
the  spacing  d/c,  the  stagger  angle  x»  and  the  Interblade  phase  angle  which 
depends  on  the  ratio  of  the  wave  length  of  the  disturbance  to  the  cascade 
spacing. 

The  essential  features  of  the  3-D  gust  interacting  with  a  linear  cascade  are: 

•  The  resonance  effect  that  is  the  conditions  at  which 
the  unsteady  lift  vanishes.  This  effect  occurs  even 
at  very  large  spacing  as  shown  In  Figure  14.  As  the 
Mach  number  Increases,  the  cascade  lift  suddenly 
collapses  while  that  of  a  single  airfoil  Is  not  affected. 

•  The  cascade  effect  Is  more  pronounced  at  large  Mach  number. 
This  effect  Is  Illustrated  In  Figure  15  showing  the  unsteady 
pressure  for  an  airfoil  and  a  cascade  at  M  ■  0.8. 
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2.  Nonlinear  Theory  for  a  2-D  Gust 

Figure  16  shows  schematically  a  lifting  airfoil  subject  to  a 
transverse  and  longitudinal  gust.  For  this  problem  the  Important  parameters 
are  the  reduced  frequency  kj,  the  transverse  wave  number  k2,  as  well  as  the 
mean  flow  round  the  airfoil.  The  most  Interesting  physical  feature  of  this 
problem  Is  the  coupling  between  the  mean  flow  and  the  oncoming  vortical 
disturbance.  This  coupling  distorts  the  wave  length  of  the  gust  as  It 
Interacts  with  the  airfoil.  The  theory  was  developed  by  Goldstein  and  Atassl 
[8],  and  more  detailed  results  for  cambered  airfoils  at  finite  incidence  to 
the  mean  flow  were  given  by  Atassl  [23,24]. 

Figure  17  shows  the  variation  of  the  real  and  imaginary  parts  of  the  lift 
function  L  for  three  different  airfoils  having  the  same  mean  loading.  The 
reduced  frequency  Is  varied  from  0  to  10,  and  the  gust  angle  was  taken  to  be 
45°  (ki  *  k2).  The  magnitude  F  of  the  ratio  of  the  unsteady  lift  acting  on 
the  airfoils  of  Figure  17,  to  the  Sears  function  Is  plotted  in  Figure  18 
versus  the  reduced  frequency.  This  clearly  shows  the  significant  nonlinear 
effects.  One  can  single  out  the  following  nonlinear  effects  on  the  gust 
response  function: 


(1)  The  coupling  between  the  mean  potential  flow  and  the 
unsteady  oncoming  vortical  disturbance  depends  In  a 
complex  manner  on  the  meanflow  Incidence,  the  geometry 
of  the  airfoil  and  the  gust  angle.  It  is  not  possible 
to  characterize  these  effects  by  a  simple  parameter  such 
as  loading.  This  Is  clearly  seen  In  Figures  17  and  18 
where  we  took  three  airfoils  having  the  same  loading. 

(11)  The  effect  of  the  longitudinal  gust  Is  strongest  at  low 
reduced  frequency.  Hie  unsteady  lift  is  significantly 
reduced  because  of  the  airfoil  loading. 

(Ill)  At  high  reduced  frequency,  the  response  function  does  not 
decay  but  seems  to  tend  toward  a  periodic  function  of 
finite  amplitude. 
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3.  Nonlinear  Theory  for  a  3-D  Incompressible  Gust 
For  an  Incompressible  flow  Eq.  (9)  reduces  to 

v2$  *  -  7*u  (13) 

and  hence  the  problem  reduces  simply  to  solving  Poisson's  equation.  Note 
however  that  because  of  vortex  shedding  In  the  wake,  the  solution  will  have  a 
discontinuity  A*  ,  whose  expression  we  have  derived  In  [2]. 

In  our  review  of  linear  theories  we  noted  that  an  Integral  formulation 
has  been  used  for  a  single  airfoil  and  a  cascade.  The  numerical  calculations 
can  be  significantly  simplified  by  transforming  the  differential  governing 
equations  Into  singular  Integral  equations.  In  fact  for  a  cascade  of 
osclllaltng  airfoils,  we  have  used  such  an  Integral  formulation  in  [25-29]. 

For  the  present  problem  Eq.  (13)  can  also  be  transformed  into  an  Integral 
equation.  Details  will  be  given  In  [30],  and  lead  to 

♦  *  "  Tfif-  ///  -  70  7G  <JV 

R 

'  *F"  ff  *  dS  +  ff  lw  dSw  *  (14) 

where  u^R)  and  4.  are  functions  determined  using  our  splitting  theorem  and 
G  Is  a  Green's  function.  The  two  surface  Integrals  In  (14)  are  over  the 
airfoil  (cascade)  surface  B  and  the  wake  sheet  W. 

Note  that  If  we  can  determine  the  Green's  function  G  for  (13)  such  that 

*  0  along  B 

then  4  Is  given  In  terms  of  surface  and  volume  integrals  of  already  known 
functions. 
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The  first  problem  and  certainly  the  easiest  to  which  the  present  theory 
Is  applied  Is  the  case  of  3-D  gust  acting  upon  an  airfoil  of  Infinite  span. 

In  this  case,  the  potential  function  $  of  the  mean  velocity  and  the 
streamf unction  t  can  be  readily  determined.  We  consider  Joukowskl  airfoils  In 
the  present  study.  This  class  of  airfoils  Is  very  popular  and  the  expressions 
for  t  and  <|>  can  be  obtained  In  closed-form  analytically. 

First  we  have  to  construct  our  regular  vortical  solution  u^RK  We  take 

4 

for  Lagranglan  coordinate  a  Llghthlll's  dri ft  function  A  ,  the 
streamf unction  ^  and  the  ordinate  z  In  the  span  direction,  a  was  defined 
in  [31]  and  Its  difference  between  two  points  on  a  streamline  Is  equal  to  the 
time  It  takes  a  fluid  particle  to  traverse  the  distance  between  those  points. 

(1)  Numerical  Scheme 

In  the  present  case,  the  Integral  equation  (14)  reduces  to  the 

form 

♦  -  ir$  *  H  (15) 

where  the  operator 

3K 

dS,  (16) 

A  dn 

4  4  4 

where  r  *  lx  -  x0|,  dS  *  |dx0|,  K0  Is  the  modified  Bessel  function  of  order 

zero,  and  Is  the  derivative  In  the  normal  function  to  the  airfoil 

surface.  The  contour  Integral  Is  taken  along  the  airfoil  A.  The 
Inhomogeneous  term  H  contains  a  double  integral  and  a  line  Integral  of  known 
functions  and  depends  on  an  arbitrary  constant. 


Two  conditions  will  be  Imposed  on  the  solution: 

•  The  solution  should  satisfy  Kelvin's  theorem. 

•  The  solution  should  satisfy  the  Kutta  condition. 

The  first  condition  can  easily  be  satisfied  by  imposing  the  requirement  that  <|> 
Is  continuous  at  the  trailing  edge.  The  second  condition  means  that  the 
pressure  Is  continuous  at  the  trailing  edge.  It  can  be  written  as 

Y  u)  =  Adk^  +  ifl0i  ||  3  -  cv  =  o  U7) 

where  Cv  Is  a  constant  depending  on  the  vortical  solution  of  the  gust  wave. 

As  we  examine  Eq.  (15),  we  see  that  the  homogeneous  equation 

X  U)  -  *  0  (18) 

has  an  eigenfunction  $0  corresponding  to  the  eigenvalue  w  of  the  operator 

X  .  and  hence  the  solution  to  (15)  Is  not  unique.  However,  the  solution 

will  be  unique  by  Imposition  of  condition  (17).  From  the  practical  point  of 

view,  a  discretization  of  (15)  will  lead  to  a  singular  system  of  linear 

algebraic  equations.  If  we  relax  one  of  these  equations  and  replace  It  by 

(17),  the  resulting  system  will  still  be  close  to  singular,  and  very  difficult 

to  solve  with  precision.  This  situation  Is  common  when  one  attempts  to  solve 

equations  where  the  desired  solution  Is  multivalued  to  account  for  the  lift. 

However,  usually  the  equation  Is  Laplace's  equation,  and  instead  of  K0  we 

have  tnr.  In  this  case  the  eigenfunction  $0  corresponding  to  it  Is  a 

constant.  In  [27],  we  were  able  to  use  this  property  to  solve  the  system  of 

linear  equations.  Here  a  more  elaborate  scheme  Is  required. 

1.  Solve  the  homoieneous  eigenvalue  for  the  Hermltian  operator. 

* 

Let  to  be  such  an  el gnef unction. 
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2.  Solve  the  modified  integral  equation. 

{at  -  ♦**>♦  a  H  (19) 

It  can  be  shown  that  (19)  Is  equivalent  to  (16)  and  (17),  and  it  has  a  unique 
solution.  This  Integral  equation  has  been  solved  by  collocation.  The  airfoil 
contour  has  been  divided  into  N  sections  and  the  value  of  the  potential 
function  has  been  assumed  to  be  constant  over  each  section.  Then  the 
requirement  for  the  Integral  equation  to  be  satisfied  in  a  set  of  N 
mid-section  points  gives  us  the  system  of  N  linear  algebraic  equations  for  N 
unknown  values  of  the  potential  function.  Also  numerical  representation  of 
the  Kutta  condition  Is  needed.  This  has  been  done  by  replacing  the  value  of 
the  derivative  of  the  potential  function  of  the  trailing  edge  by  a  ratio  of 
two  finite  differences.  Such  approximation  requires  closer  spacing  of  the 
collocation  points  at  the  trailing  edge  region.  Similar  spacing  has  been 
applied  also  for  the  leading  edge  region  where  potential  function  gradients 
are  unusually  large. 

The  unsteady  pressure  along  the  airfoil  surface  is  given  by 

i  i  D 

P  “  PTE  *  '  p  fit  (*} 

This  expression  contains  the  derivative  of  the  potential  function.  However, 
the  potential  function  Is  known  only  in  a  discrete  set  of  points  of  the 
airfoil  surface.  Finite  difference  behavior  for  the  derivative  approximation 
cannot  be  used,  because  It  leads  to  discontinuous  pressure  distribution. 
Therefore  spline  Interpolation  of  the  potential  function  has  been  made  which 
gives  continuous  derivatives  along  the  airfoil  surface.  The  first 
calculations  showed  that  for  the  third-order  polynomial  spline  much  more 
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collocation  points  are  needed  to  insure  considerable  smoothness  of  the 
pressure  distribution  for  the  desired  accuracy.  For  that  reason  the 
exponential  spline  has  been  used  instead  of  the  polynomical  spline. 

(ii)  Results 

As  an  application  of  our  theory  we  have  chosen  a  typical  Joukowski 
airfoil  with  camber,  thickness,  and  incidence  to  the  mean  flow.  The  airfoil 
geometry  is  shown  in  Figure  4.  It  has  0.1  camber,  0.1  thickness,  and  is 
placed  at  10°  incidence  to  the  mean  flow.  This  airfoil  is  subject  to  a 
sinusoidal  gust  characterized  by  three  wave  numbers,  the  reduced  frequency, 
ki  =  1.0,  k2  =  1.0,  and  k3  =  1.0 

Drift  Function 

The  drift  function  4  was  used  as  a  Lagrangian  coordinate  along  a 
streamline.  This  enabled  us  to  integrate  the  vorticlty  equation.  However,  a 
which  represents  the  time  increment  as  we  move  along  a  streamline,  is  rather  a 
complicated  function  of  the  fluid  particle  position.  Figure  19  shows  a 
three-dimensional  plot  of  0  =  a  -  <(»0  ,  where  4>0  is  the  mean  potential  velocity 
function  versus  the  coordinate  x  and  the  streamf unction  i|>  .  The  lines 
stretching  from  left  to  right  correspond  to  constant  x,  while  those  stretching 
in  the  x-direction  correspond  to  constant  (along  the  same  mean  streamline). 
Note  that  far  from  the  airfoil,  located  between  -1  and  +1,  there  are  two 
plateaux  for  D  at  different  levels  because  of  the  circulation  round  the 
airfoil.  At  streamlllnes  closer  to  the  suction  side  of  the  airfoil  there  is  a 
significant  decrease  In  the  value  of  D  until  we  get  very  close  to  the  airfoil, 
then  the  effect  of  the  stagnation  point  begins  to  be  felt  and  as  a  result  D 
Increases  sharply  to  become  Infinite  at  the  stagnation  point  and  remains 
Infinite  along  the  airfoil  surface  and  its  wake.  On  the  other  hand  at 
streamlines  close  to  the  pressure  side  of  the  airfoil  surface  D  Increases 


monotonically  as  we  move  along  a  streamline.  It  also  Increases  as  we  move 
closer  to  the  airfoil  until  it  becomes  Infinite  at  the  stagnation  point  and 
the  airfoil  surface  and  its  wake.  The  case  <i>  *  0  corresponding  to  the  airfoil 
surface  is  not  plotted  because  of  the  singular  behavior  of  D. 

It  is  this  behavior  of  the  drift  function  that  produces  the  singular 
behavior  of  the  vortlcity  at  the  airfoil  surface.  The  regular  vortical 
velocity,  however,  does  not  have  this  singular  behavior. 

The  Vortical  Velocity  u(R) 

The  expression  of  the  vortical  velocity  u^R)  is  determined 
analytically.  The  analytical  expression,  however,  depends  on  the  drift 
function  and  the  potential  and  stream  function  $0  and  \|>0  .  These  three 
quantities  are  explicitly  known  only  numerically  in  terms  of  Cartesian 
coordinates.  In  order  to  Illustrate  the  result  of  our  regularization 
procedure  we  have  plotted  in  Figures  20  and  21  the  magnitudes  of  the 
streamwlse  and  normal  components  of  u^R),  versus  x  and  i|>0  respectively. 

Again,  the  lines  stretching  from  left  to  right  correspond  to  constant  k,  while 
those  stretching  in  the  x-directlon  correspond  to  constant  (along  the  same 
mean  streamline).  These  figures  clearly  show  that  both  the  streamwlse  and 
normal  components  of  u^R)  vanish  at  the  airfoil  surface  as  predicted  by  the 
theory.  The  normal  component  Un  exhibits,  however,  very  strong  gradient  near 
the  airfoil  surface,  while  the  streamwlse  component  has  a  rather  smooth 
behavior. 

The  Potential  Function  $ 

The  function  4  is  a  solution  to  the  integral  equation  (14)  which 


after  modification  became  equation  (19).  Figures  22  and  23  show  plots  of  the 
real  and  Imaginary  part  of  $  versus  the  airfoil  arc  length  counted  from  the 
pressure  side  trailing  edge.  The  vertical  solid  lines  at  S  *  1.97  and  2.02 
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correspond  to  the  locations  of  the  stagnation  points.  Note  that  <j>  is 
essentially  a  smooth  curve,  with  strong  gradients  only  in  the  stagnation  point 
area. 

The  Unsteady  Pressure 

The  unsteady  pressure  is  the  material  derivative  of  $.  Figures  24 
and  25  show  plots  for  the  real  and  imaginary  components  of  the  unsteady 
pressure  for  two  different  airfoils  also  plotted  on  these  figures.  For  a 
thin  airfoil  at  zero-incidence  (Fig.  24),  the  unsteady  pressure  is  most 
important  near  the  leading  edge.  This  result  is  similar  to  that  of  the  flat 

plate  case  where  the  pressure  is  infinite  at  the  leading  edge.  For  a  10% 

thickness  airfoil,  the  real  component  of  the  pressure  in-phase  with  the  gust 
velocity  does  not  exhibit  large  values  near  the  leading  edge,  while  the 
imaginary  out-of-phase  component  has  a  large  peak  in  this  area. 

It  is  of  course  not  possible  to  draw  a  conclusion  about  the  behavior  of 
the  pressure  from  the  present  results.  More  studies  are  needed  to  determine 

the  essential  features  of  the  unsteady  pressure  in  terms  of  the  many 

parameters  entering  the  problem. 

The  Unsteady  Lift 

The  present  theory  is  the  most  general  aerodynamic  analysis  for  an 
arbitrary  shaped  airfoil  subject  to  three-dimensional  gust.  It  contains  all 
previously  studied  cases  such  as  that  of  Sears  [3],  the  nonlinear  2-D  second 
order  theory  [8],  and  Graham  [22]  oblique  gust  theory  for  a  flat  plate 
airfoil.  Therefore  we  started  by  comparing  our  results  with  the  linear  theory 
of  Graham  and  our  second  order  theory. 

We  first  calculated  the  real  and  imaginary  parts  of  the  lift  for  a  flat 
plate  airfoil  subject  to  an  oblique  gust  with  a  spanwlse  wave  number  k3  *  1.0. 
Figure  26  shows  the  imaginary  part  of  the  lift  (out-of-phase  component)  versus 
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the  real  part  of  the  lift  (in-phase  component)  for  different  values  of  the 
reduced  frequency  Iq.  Because  we  cannot  take  a  zero- thickness  airfoil,  we 
have  considered  a  $%  thickness  airfoil.  Our  results  are  plotted  in  a  similar 
way  to  the  flat  plate  case  and  indicate  an  excellent  agreement  up  to  reduced 
frequencies  of  order  2.  Above  kj  =  2,  a  slight  difference  in  the  phase  of  the 
lift  appears  to  exist  between  the  linear  theory  and  ours.  This  could  be 
attributed  to  the  5%  thickness  of  our  airfoil.  However,  it  could  also  result 
from  a  reduced  accuracy  in  the  numerical  scheme  as  kj  increases,  to  larger 
values.  A  further  investigation  of  the  high  frequency  case  is  underway  to 
clarify  this  point. 

The  second  comparison  we  have  carried  out  was  with  our  previous  results 
[24],  In  [24]  we  have  derived  an  analytical  formula  for  the  lift  of  a 
cambered  airfoil  at  incidence  to  the  mean  flow  and  subject  to  a 
two-dimensional  gust.  Figure  27  shows  the  result  of  the  second  order  theory 
(SOT)  compared  with  our  present  (PRT)  calculation  for  a  5 %  thickness  airfoil. 
Again,  we  see  an  excellent  agreement  up  to  a  reduced  frequency  of  order  2. 
Figure  28  shows  similar  results  for  the  same  airfoil  but  with  5°  angle  of 
attack . 

These  comparisons  show  clearly  the  good  accuracy  of  our  computed  results. 
Figure  29  shows  the  three-dimensional  effect  of  the  gust,  where  we  have 
compared  the  lift  resulting  from  a  2-D  gust  with  that  resulting  from  a 
three-dimensional  gust  k3  =  1,  the  reduced  frequency  kj  was  varied  from  0.25 
to  2.5  .  The  gust  transverse  wave  number  is  k2  ■  V.\  .  It  is  noted  that  k3 
has  a  very  significant  effect  on  both  the  magnitude  and  the  phase  of  the  lift 
low  frequency.  In  fact  for  ki  between  0.25  and  0.5,  both  magnitude  and  phase 
of  the  3-0  gust  are  almost  constant.  This  result  is  In  harmony  with  certain 
results  given  by  Graham  [25]  and  in  fact  should  be  expected  in  view  of  the 
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fact  that  the  airfoils  considered  departs  only  very  slightly  from  a  flat 
plate  airfoil. 

The  comparisons  we  have  carried  out  with  the  2-D  linear  Sears  results, 

3-D  linear  Graham  calculations,  and  our  previously  developed  2-D  nonlinear 
theory,  are  not  yet  complete.  For  example,  it  Is  very  Important  to  determine 
the  validity  of  the  2-D  nonlinear  theory  as  the  mean  lift  of  the  airfoil  is 
Increased.  This,  of  course,  can  now  be  done  by  a  more  detailed  comparison  of 
our  present  work  which  encompasses  all  previous  cases  as  just  a. particular 
case  assigned  to  certain  parameters.  An  Investigation  comparing  all  these 
theories  is  being  carried  out. 

V.  COMPARISON  BETWEEN  THEORY  AND  DATA 
FOR  HIGHLY  LOADED  OSCILLATING  CASCADES 

We  have  previously  developed  under  AFOSR  sponsorship  an  aerodynamic 
theory  for  oscillating  airfoils  In  cascade  with  arbitrary  airfoil  and  cascade 
geometry.  The  airfoils  could  have  translational  and  rotational  oscillations 
with  circumferential  modes  modeled  as  an  Interblade  phase  angle  In  the 
oscillatory  motion.  The  theory  was  published  in  [27].  The  analysis  was  later 
applied  to  calculate  turbine  blade  flutter  In  [28]  and  [29],  The  results 
clearly  show  that  coupled  bending  and  torsional  modes  could  produce  flutter  at 
typical  operating  conditions  for  highly  loaded  blades.  In  an  "Overview  of 
NASA/AF/NAVY  Symposium  on  Aeroelastlclty  of  Turbine  Engines  [32],  Professor 
Slsto  noted  the  importance  of  this  result  and  stated  that  "Experimental 
confirmation  of  the  validity  of  this  conclusion  should  be  anticipated  by  the 
R&D  community,"  [p .7] . 

At  UTRC,  Dr.  Franklin  Carta  has  been  Investigating  the  pressure 
distribution  on  oscillating  cascades  In  pitching  oscillation.  Dr.  Carta's 


research  is  sponsored  by  NASA  Lewis  Research  Center.  The  above  results 
encouraged  us  to  cooperate  for  a  comparison  between  theory  and  experiment. 

Or.  Carta  visited  us  at  Notre  Game  and  provided  us  with  details  regarding  his 
experimental  condition.  The  geometry  of  his  airfoils  and  the  cascade 
parameters  are  given  in  Table  1.  We  carried  out  a  study  of  his  cascade  for 
interblade  phase  angles  varying  from  0  to  360®.  Similar  comparisons  were  also 
carried  out  with  Verdon  and  Caspar  [33],  also  from  UTRC,  using  a  purely 
numerical  code.  The  results  were  presented  in  [34]. 

Figure  [30]  shows  essentially  this  comparison  for  a  reduced  frequency, 
k  a  0.122  and  four  values  of  the  interblade  phase  angle  a  .  The  results  show 
excellent  agreement  between  our  analysis,  Verdon  and  Caspar  numerical  code, 
and  the  experimental  data.  Note,  however,  the  strong  deviation  of  the  real 
part  of  the  pressure  in  the  Verdon/Caspar  results  near  the  trailing  edge. 

This  deviation  is  caused  by  the  difficulty  of  accurately  capturing  the 
singular  behavior  in  the  unsteady  pressure  at  a  sharp  trailing  edge  with  a 
finite  difference  scheme. 

VI.  FUTURE  PLANS 

Our  past  and  present  work  under  AFOSR  sponsorship  has  always  covered  a 
broad  range  of  topics  which  have  evolved  from  the  initially  stated  objectives. 
This  is  particularly  true  for  the  period  covered  by  the  present  report. 

(i)  We  started  out  to  study  the  aerodynamics  of  airfoils  subject  to 
three-dimensional  unsteady  vortical  disturbances.  A  problem  of 
great  interest  in  aerodynamics  and  aeroelastlclty. 

(11)  This  led  us  to  develop  new  mathematical  methods  which  apply  to  large 
structure  disturbances  acting  upon  streamlined  bodies  with  wakes 
such  as  airfoils  and  blades  in  turbomachines  subject  to  rotor-stator 
Interaction  or  inlet  distortion. 

(Ill)  Later  we  generalized  the  mathematical  procedure  to  cover  the  equally 
Important  cases  of  small  structure  turbulence  disturbances 
interacting  with  streamlined  and  bluff  bodies.  This  is  known  as  the 
rapid  distortion  theory  of  turbulence. 
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( 1 v)  The  mathematical  theory  led  to  a  uniformly  valid  splitting  theorem 
for  unsteady  vortical  and  entropic  disturbances  of  potential  flows. 

(v)  This  led  us  to  Initiate  an  Investigation  to  study  the  change  in 

turbulance  near  the  stagnation  point  of  a  flow.  A  problem  of  great 
Interest  in  studying  the  stability  of  a  flow  as  It  Interacts  with  a 
body. 

The  above  shows  the  broad  scope  of  our  research  interests  which  essentially 
are  In  the  following  areas: 

•  Unsteady  aerodynamics  of  airfoils  and  cascades 

•  Flutter  and  stability  analysis  of  airfoils,  cascades 

and  assembled  systems 

•  Mathematical  methods:  analytical  and  numerical 

•  Unsteady  flows 

Our  future  plans  are  therefore  determined  by  ongoing  state-of-the-art 
developments  in  the  aerodynamics  and  the  stability  analysis  of  turbomachine 
systems  and  their  relevance  to  technological  developments. 

Essentially,  these  are  the  problems  we  shall  be  pursuing  and  developing 
in  the  near  future: 


1.  Three-Dimensional  Vortical  Disturbances  Acting  Upon  Loaded  Airfoils 

We  shall  complete  the  analytical  analysis  of  the  aerodynamics  of  an 
airfoil  of  arbitrary  geometry  subject  to  3-D  gust.  This  will  include  a 
parametric  study  of  the  unsteady  forces  In  terms  of  the  airfoil  geometry, 
the  mean  flow  Incidence,  and  the  frequency  and  the  direction  of  the  gust 
as  defined  by  the  three  wave  numbers.  The  comparison  with  linear  3-D 
theories  and  nonlinear  2-D  theories  will  also  be  completed.  This  study 
should  determine  the  significant  parameters  affecting  the  aerodynamics  of 
a  3-D  gust. 

2.  Interaction  of  Swirling  Vortices  with  Loaded  Airfoils 

In  turbomachines,  major  upstream  disturbances  are  produced  either  by 
viscous  wake  defects  or  by  swirling  vortices  usually  emanating  from 
upstage  tips.  It  Is  Important  for  the  designer  to  know  which  effect 
produces  higher  fluctuating  lift  and  whether  a  design  shape  optimization 
can  reduce  the  unwanted  fluctuations  In  the  forces.  This  study  can  be 
carried  out  as  an  application  to  our  general  theory  of  a  3-D  gust. 
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3.  High  Frequency  3-D  Disturbances  Acting  Upon  Loaded  Airfoils  at 
High  Subsonic  Mach  Numbers- 

This  Is  the  case  closest  to  operating  conditions  for  high-speed 
machines.  The  treatement  of  the  general  problem  requires  significant 
mathematical  and  numerical  efforts.  However,  If  one  assumes  the  reduced 
frequency  to  be  high,  a  simplification  of  the  mathematical  equation 
governing  the  problem  can  be  brought  about  and  leads  to  a  formulation 
similar  to  the  one  we  have  developed  for  a  low-Mach  number. 

4.  3-D  Vortical  Disturbances  Acting  Upon  Loaded  Cascades 

Depending  on  the  solidity,  the  cascade  effects  can  be  very 
Important.  Even  for  a  small  solidity,  resonance  effects  can  take  place 
In  a  cascade.  The  extension  of  our  theory  to  a  cascade  of.  airfoils  can 
be  done  with  no  difficulty  except  for  the  added  computational  time  due  to 
define  the  mean  flow  of  the  cascade.  The  incorporation  of  high  Mach 
number  and  high  frequency  to  the  cascade  analysis  will  be  a  very 
Important  contribution  to  a  more  accurate  evaluation  of  flutter  and  noise 
generation  in  turbomachines. 

5.  Change  in  Turbulence  and  Pressure  Fluctuation  Near  a  Stagnation  Point 

As  a  spin-off  of  our  splitting  theorem,  we  can  study  the  change  in 
the  characteristics  of  a  turbulent  flow  near  the  stagnation  point  of  an 
airfoil  or  a  bluff  body.  This  Is  Important  to  determine  the  stability  of 
a  flow  as  It  interacts  with  the  body.  Separation  and  boundary- layer 
stability  are  affected  by  the  level  of  turbulence  and  the  pressure 
fluctuation  near  a  stagnation  point. 

6.  Aeroelastic  Characteristcs  of  Loaded  Blade  Disc  Assemblies 


All  present  aeroelastic  stability  studies  of  blade  disc  assemblies 
use  simplified  aerodynamic  codes  developed  from  linear  unsteady 
aerodynamic  theories  wherein  the  blades  are  approximated  by  flat  plates. 
Our  previously  developed  theory  for  loaded  cascades  oscillating  with 
circumferential  modes  showed  significant  changes  in  the  stability 
boundary  from  the  results  of  flat  plate  cascades.  We  feel  that  a  study 
Is  needed  to  Incorporate  a  more  realistic  aeroelastic  model  with  our 
aerodynamic  theory.  Such  a  study  will  define  a  lore  accurate  flutter 
boundary  for  loaded  blade  disc  assemblies. 
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Angle  of  Attack  a  =  0.0785398  2  4.5° 

THICKNESS  =  0.0601266  Chord 

CAMBER  =  0.0218513  Chord 

GAP/CHORD  =  0.7500880 

STAGGER  =  0.959930 

REDUCED  FREQUENCY  =  0.122000 
o  =  INTERBLADE  PHASE  ANGLE 

TABLE  1.  United  Technologies  Cascade  Parameters 


FIGURE  1.  SCHEMATIC  REPRESENTATION  OF  FLOW  BEHIND  A  ROTOR  BLADE  ROW 


ROTOR  WAKES  INTERACTION  WITH 
DOWNSTREAM  STATOR.. 


FIGURE  3.  ROTOR  WAKES  INTERACTION  WITH  DOWNSTREAM  STATOR 


STREAMWISE  COMPONENT  OF  THE  VORTICITY  (IMAGINARY) 
JOUKOWSKI  AIRFOIL:  CAMBER=0 . 1  THICKNESS=0 . 1 


A2-0.71 


FIGURE  8.  THE  NORMAL  COMPONENT  OF  THE  VORTICITY  (IMAGINARY) 
JOUKOWSKI  AIRFOIL:  CAMBER=0 . i  THICKNESS=0 . 1 
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FIGURE  13.  REAL  AND  IMAGINARY  PARTS  OF  THE  LIFT  COEFFICIENT  L  FOR  AN  OBLIQUE 
GUST  ACTING  ON  A  FLAT  PLATE  FOR  DIFFERENT  SPANWISE  WAVE  NUMBER  k3 
THE  REDUCED  FREQUENCY  k,  IS  VARIED  FROM  0  to  5. 
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MACH  NUMBER 


FIGURE  14.  COMPARISON  BETWEEN  THE  LIFT  OF  A  SINGLE  AIRFOIL 
AND  A  CASCADE  OF  VERY  LARGE  SPACING.  RESONANCE 
OCCURS  AT  M  =  0.931. 
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FIGURE  15.  COMPARISON  BETWEEN  THE  AERODYNAMIC  PRESSURE 
DISTRIBUTION  ALONG  A  SINGLE  AIRFOIL  AND  A 
CASCADE  BLADE.  THE  INTERBLADE  PHASE  ANGLE  =  2. 


REAL  COMPONENT  OF  LIFT 


OF  THE  UNSTEADY  LIFT  FOR  LIFTING  AIRFOILS  TO  THE 
REDUCED  FREQUENCY. 


FIGURE  19.  THE  DRIFT  FUNCTION  MINUS  POTENTIAL  FUNCTION 
JOUKOWSKI  AIRFOIL  (CAMBER-O.l  THICKNESS-0 . 1  INCIDENCE-10  DEG) 


.71 


FIGURE  20.  THE  MAGNITUDE  OF  THE  STREAMWISE  VELOCITY  (HA-METHOD) 

JOUKOWSKI  AIRFOIL  (CAMBER-0 . 1  THICKNESS-O . 1) 


POTENTIAL  FUNCTION:  REAL  (SOLID  LINE)  AND  IMAGINARY  (DASHED  LINE) 
JOUKOWSKI  AIRFOIL:  CAMBER-0 . 05  THICKNESS-0.05 


Al  — 0.71  Kl-i 
A2-0.71  K2-1 

A3-0  K3-1 

INCIDENCE-10  DEG 


POTENTIAL  FUNCTION:  REAL  (SOLID  LINE)  AND  IMAGINARY  (DASHED  LINE) 
JOUKOWSKI  AIRFOIL:  CAMBER-0.1  THICKNESS-0.1 
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Al=-0 .71  K 1= 1 
A2-0.71  K2=l 
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FIGURE  24.  UNSTEADY  PRESSURE 

REAL  (SOLID  LINE)  AND  IMAGINARY  (DASHED  LINE) 

J0UK0WSKI  AIRFOIL:  CAMBER-0.05  THICKNESS-0.05 
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FIGURE  26.  THE  UNSTEADY  LIFT  COEFFICIENT 

JOUKOWSKI  AIRFOIL:  CAMBER=0  THICKNESS=0 . 05 
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E  28.  THE  UNSTEADY  LIFT  COEFFICIENT 
JOUKOWSKI  AIRFOIL:  CAMBER=0  THICKNESS=0 . 05 


JRE  29.  THE  UNSTEADY  LIFT  COEFFICIENT 
JOUKOWSKI  AIRFOIL:  CAMBER=0  THICKNESS 
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FIGURE  30.  COMPARISON  OF  TWO  THEORIES  WITH  EXPERIMENT  FOR 


VIII.  PUBLICATIONS  RESULTING  FROM  OUR  WORK  UNDER  AFOSR  SPONSORSHIP 


All  publications  listed  below  are  authored  by  Professor  H.  Atassi.  Only 
co-authors  are  listed. 


"Three-Dimensional  Periodic  Disturbances  Acting  Upon  Airfoils  in  Cascade", 
Aeroelasticity  in  Turbomachines,  Ed.  P.  Suter,  Juris-Verlag  Zurich,  pp. 
383=398  ,"1981" - - - 

"Stability  and  Flutter  Analysis  of  Turbine  Blades  at  Low  Speed", 
Aeroelasticity  in  Turbomachines,  Ed.  P.  Suter,  Juris-Verlag  Zurich,  pp. 
187-201 , 1981.  Co-author  T.J. Akai . 


In  Preparation 

"Regularization  of  Goldstein's  Splitting  of  Unsteady  Vortical  and  Entropic 
Distortions  of  Potential  Flows,"  Invited  paper,  19th  Annual  Meeting,  Society 
of  Engineering  Science,  October  27-29,  1982,  Rolla,  Missouri. 

"A  Uniformly  Valid  Splitting  of  Unsteady  Vortical  and  Entropic  Disturbances  of 
Potential  Flows." 

"Three-Dimensional  Periodic  Vortical  Distrubances  Acting  Upon  an  Airfoil," 
Co-author  J.  Grzedzinski 


IX.  INVITED  LECTURES 


These  are  lectures  given  by  Professor  H.  Atassi  on  topics  related  to  the 
present  AFOSR  grant. 

"Unsteady  Aerodynamics  of  Lifting  Airfoils,"  NASA  Lewis  Research  Center, 
Cleveland,  Ohio,  December  10,  1981. 

"Unsteady  Aerodynamics:  Early  and  Recent  Developments,"  Ohio  State  University, 
Columbus,  Ohio,  May  21,  1982. 
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